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1, INTRODUCTION 

Public transportation nowadays is one of the most encountered daily problems in most cities of the 
world. Due to many factors including world fuel crises, environmental concerns, there was a demand for fuel 
economy and environment friendly transportation vehicles [1]. Electric Vehicles (EV) and Electric Transport 
Systems (ETS) had acquired great attention and are predicted to be the most used mode of transportation in 
the near future, due to their advantages as compared to ordinary transportation vehicles [2]. These vehicles 
are economic with low running coasts and maintenance, have better performance capabilities as compared to 
conventional internal combustion engine vehicles. Since these vehicles use no fuel, they are predicted to 
reduce harmful carbon emissions by 25-60% and thus reducing pollution and greenhouse gases. Moreover, 
the use of renewable energy for charging is an important option that helps in reducing fuel problems. In 
comparison with an internal combustion engine, electric motor has high efficiency and lower energy losses. 
The efficiency of the electric motor is about 90% against 25% of an internal combustion engine. Researchers 
are continuously working on ETS and EV to improve their range, performance and reliability [3]. In this 
field, the ongoing research efforts are focused on increasing the life time of the battery bank [4] that provides 
the advantages of low power rated converter and better load profile for the battery bank. Clegg (1996) [5] 
indicated that the energy efficiency can be amplified by using regenerative braking of motors (tractive 
motor) [6]. His results indicated a good dynamic performance and robust stability with more recovered 
energy. Accordingly, regenerative braking appears to be suitable for hybrid electric vehicles and for fuel cell 
electric vehicles as well [7]. 
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Recently, in the Capital city of Jordan, Amman, a two lane Express Bus Project of 25 km long has 
been lunched [8]. The express bus project, in its current form Figure 1 employs the use of fossil energy 
sources (Diesel). 





Figure 1. A view of Amman Express Bus Project 


However, the current investigation suggests the use of electric buses and the better investment of the 
French funding [9] in establishing the electrical network required for this project. Moreover, the study also 
suggests the use of traction electric motors instead of internal combustion engine motors that requirs fossile 
fuel (diesel), taking into consideration the geographical nature of Amman and the possibility of employing 
solar energy to generate the electricity required for running these electric buses. 
The main requirements of a traction electric motor are: 
a) The ability to realize the maximum torque for the entire speed range. 
b) The ability to work in two directions without any additional devices. 
c) The simplicity of the air cooling design. 
d) The ability to operate in the generator mode. 
Electric vehicles use several electric motors to drive individual wheels in order to increases the 
traction power of the vehicle. A series excited DC electric motors are used as traction motors [10]. 


2. THE MATHEMATICAL MODEL 

Computer modeling has always been used for industrial designs, especially for EVs which have 
benefited greatly from the accelerated use of simulation programs in evaluating all vehicles aspects, 
including energy consumption. During the current study, a mathematical model of a direct current electric 
(DC) motor with series excitation - traction electric motor - for dynamic modes, Figure 2 was established. 
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DC Supply 


Figure 2. A schematic diagram of a direct current electric motor with series excitation 


The implementation of the second Kirchhoffs’ law and the basic equation of motion give the 
possibility of obtaining the mathematical model for a DC Motor [11]. 
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Where; = K@q is the internal induced voltage of the armature. The total resistance of the motor (Rio:) 


could be calculated from the equation: 


Riot ae I pai +r IN dhes a I teppales 
The total inductance of the motor (Lio) can be easily calculated from the equation: 
Lot = ie a Levies as La arpoies 


Substituting the expression of the internal induced voltage of the motor into the second Kirchhoff 
law, gives: 


7 
i Rior + Lioy oe + Kdo =V, (2) 
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Then, supplementing the mathematical model of the DC motor with the basic equation of motion of the 
electric drive can be illustrated with the following equations: 


T-T, = (3) 
dt 


T =KGi, (4) 


The substituting of expression (4) into equation (3) gives a detailed equation of motion of the 
electric drive. 


- _p _ 7 de (5) 
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Then, rewrite equations (1) and (4) to form 








“ = - (V, — 1g Riot - Ko) (6) 
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As a result of such transformations, a nonlinear mathematical model of the direct current electric 
motor with series excitation was obtained. The non-linearity of the model is related to the nonlinear 
characteristics of the machine’s magnetic circuit, which could be attributed to the nonlinear relationship 
between the flux of the machine and the field current. The obtained equations in the employed model are 
differential equations of the first order. To solve the last two equations, Euler’s’ method [12] was employed. 
Thus, equations (6) and (7) were written as the algebraic equations (8 & 9), and were then easily 
implemented with the computer application. 
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The algorithm for the solution was: 

At | (10) 
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Where Az, is the increment of the armature current at “m” step of integration; "m" is the step of 
integration; V,, is the voltage at “m” step of integration applied to the armature circuit; Ar — is the 


integration step; 1,,,, is the value of the armature current at "m-1" which refers to the previous step of 


integration; ~__, 1s the value of the angular velocity of the DC series motor on "m -1" step of integration; 
Lam a amt a Aiam (1 1) 


Where ;, 1s the increment of the armature current at "m" integration step. 


A@ 2/EE 


m J (Kiam —T hn) (12) 


Where Aq@,, is the increment of the angular velocity at the "m" integration step and 7,,, 1s the value of the 


" " 


load moment of the drive at the "m" integration step. 


O,,=0,_,+Aa,, (13) 


Equations (10) and (12) were used for the determination of the increments in the current AL in and 
the velocity AQ, . Equations (11) and (13) on the other hand, were used to calculate the current bd and the 


velocity @,, . 


3. SIMULATION RESULTS AND ANALYSIS 

In the current investigation, two simulation models were employed and their results were evaluated. 
These included Excel and Matlab/Simulink [13] moldes. The main similarity between these method of 
simulation (Excel and Matlab/Simulink) is basically the mathematical model [14]. Excel helps in 
investigating the statistical and dynamic modes of energy [15, 16]. Briefly, for simulating an electric DC 
motor with series excitation in a dynamic mode [17], values of the independent variable (tf) are introduced to 
column I of the Excel sheet, starting with the value t= 0. 


The value of the integration step - the increment - ( Af ) is determined using the empirical formula: 





At < i (14) 


The minimal time constant (7,,;, ) 1s determined by comparing the values ;,, = . Ao ,and; = Prot, 
Cc 


a 





at 
Values of L, , obtained from equation (11), are introduced to second column (B). Values of ( @,, ) as 
calculated from equation (13) are then introduced to the third column (C) of the spreadsheet. The fourth and 
fifth columns contained the calculated values of AL in and AO, obtained from equations (10) and (12), 


respectively. Excel results are better than Matlab/Simulink for investigating the statistical modes of energy. 
All calculations using the software package “Matlab” are made with variables, as with matrices [18]. 
The transient process time (t) from 0 to the end of the transient process 7,, =5z,,, 1S given in the form of a 


matrix ¢— lo - At: Tol; where, 


Ke Git (15) 
10 


Performing the process of modeling of the transient process for an electrical circuit with an 
induction coil when turned on to a DC voltage, 1s demonstrated by the following example. The used 
mathematical model in this example is a classical solution for a linear inhomogeneous differential equation of 
the first order. The parameters of the electrical circuit: 
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L=135V =20;R=105 t4,,, = L/R 3 delta = Tym, (10 5T,,, =5t :t=|0: delta: T,, |:i = Ma 


Rise Fam) 


arm. 


plot(t, 7). 


Figure 3 shows the dependence of the current in the electrical circuit with an inductance coil, when 
switched to a constant voltage in the transient process. 
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Figure 3. The dependence of the current in an electric circuit with an inductance coil 


The most effective tool for mathematical modeling of electro technical objects in dynamic modes is 
the Matlab/Simulink [10]. Figure 4 represents a structural scheme of the mathematical model of an electric 
DC motor with series excitation based on equations (5) and (6). 





Figure 4. The mathematical model of an electric DC motor with series excitation in the form of a structural 
scheme. 


Analysis of the scheme revealed a nonlinear model. The model contains three non-linear links, using 
which the motor torque (T) and the motor back EMF were determined [19]. The model of series excitated DC 
motor with a direct and soft-linear rising voltage-startup in the program of Matlab/Simulink is shown in 
Figure 5, the Matlab/Simulink Model blocks description is listed in Table 2. The motor parameters used in 
the current simulation study are listed in Table | at: A) no-load, B) 50% rated load and C) rated load. The 
employed Matlab/Simuling blocks description 1s listed in Table 2. 


Tablel. Motor Parameters 


Quantity Symbol 
Armature Resistance Ra 0.05 O 
Electromechanical Time Ta O.lc 
Constant 
Moment of Inertia J 10Kg — m? 
Load Torque Tl 200 Nm 
Armature (Supply) Va 220 V 
Voltage 


Simulation and analysis of the possibilities of traction electric motor (Ayman Y. Al-Rawashdeh) 


0 


ISSN: 2502-4752 


Table 2. Matlab/Simulink Model blocks description 


Scope Display signals- Current, Torque and Speed- generated during simulation 
Gain Multiply input by constant 

Math Function Perform mathematical function 

Product Multiply and divide scalars and nonscalars or multiply and invert matrices 
Sum Add or subtract inputs 

Constant Generate constant value 

Step Generate step function as input voltage and load 

Ramp Generate constantly increasing input voltage 

Transfer Function Model linear system by transfer function for motor circuit 

Saturation Limit input signal to the upper and lower saturation values of input voltage 
Selector Select input elements from multidimensional signal to display 











Figure 5. Matlab/Simulink mathematical model of electric DC motor with series excitation 


Figure 6 shows the dependence of the angular velocity, armature current and the moment of electric 
DC motor with series excitation during a direct start-up obtained by the simulation in the Matlab/Simulink 
program. Figure 7 shows the dependence of the angular velocity, armature current and the moment of electric 
DC motor with series excitation during a direct start-up, obtained by simulation in the Matlab/Simulink 
program. 





Figure 6. The direct start-up dependence of the angular velocity, the armature current and the moment of 
electric DC motor with series excitation at A: at no-load; B: at 50% of rated load and C: at rated load 


As could be concluded from figures 5 and 6, the simulation results revealed that the starting current 
of the motor was 460 A, but after a transient condition period, the current reached a value of 357 A and 
became stable. Similarly, at starting conditions, the motor torque reached a value of 645 N.m but after 12 s, 
the motor torque was stabilized to a value of 500 N.m [6]. The figures also indicated that the speed was 
changing smoothly (0-145) rad/s and was stable. The above mentioned values of current, torque and speed 
are for the smooth starting method [10]. 
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The energy efficiency and the minimizing of the electric energy consumption in the ETS can be 
amplified by using regenerative braking. Depending on the nature of the roads, there is a good possibility for 
the utilization of regenerative braking with good dynamic performance, more recovered energy and with a 
chance of choosing a better control algorithm for the vehicle controller [20]. When driving on a flat road 
surface at a constant speed, the consumption of electric energy is the smallest. An increase in the device 
weight, road resistance and resistances while driving up or down the hill will decrease or increase the ability 
of acceleration and thus increaseing or decreasing the energy consumption during the acceleration or de- 
acceleration. However, more energy is recovered during engine braking. Most energy is consumed during the 
acceleration of the electric transport device [21, 22]. 


A 





Figure 7. The start-up with a linear rising voltage (pulse-width modulation) dependence of the 
angular velocity, armature current and the moment of electric DC motor with series excitation: A: at no-load; 
B: at 50% of rated load and C: at rated load 


The effect of sudden changes in motor load, road resistance, passengers load or road slope were also 
investigated, as shown in Figure 8. Results revealed that the speed, torque and current were not very sensitive 
to such changes indicating that the motor performance and stability were very good. The results of the current 
simulation study indicated that the series DC motor reached a steady state speed after a transient condition 
period not exceeding 1.5 — 2 s and in different loading cases of the electric transport device. The findings of 
the current study confirms that the Series DC Motor is covering the electric transport requirements with best 
dynamic performance and is suitable to be used in electric transport devices and generally for electric 
traction. 

Figure (8) shows the dependence of the angular velocity, armature current and the moment of 
electric DC motor with series excitation during direct start-up and steady state conditions, obtained by 
simulation in the Matlab/Simulink program when the load of motor was changed during motor running, due 
to the road resistance and resistances while driving up or down the hill. 


- 


Figure 8. The direct start-up dependence of the angular velocity, the armature current and the moment of 
electric DC motor with series excitation when the load changing during motor running A: with load 
increasing (up the hill) and B: with load decreasing (down the hill) 


4. CONCLUSION 

The obtained characteristics by Simulink/Matlab of the investigated Series DC motor in terms of 
speed, current and torque are compatible with the findings of other researchers [23]. In the current work, a 
Matlab/Simulink mathematical model derivation for the described motor has been carefully described and 
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discussed taking into consideration the nonlinear relationship between the field current and the machine flux. 
Simulation results obtained in the current study confirmed the feasibility of the use of the series excited DC 
motor (traction motor) in electric buses in transportation system as alternatives for internal combustion 
engine motors. The obtained characteristics for speed, current and torque indicated that the studied motor 
reached a steady state speed after 6 s at different loads. This result clearly revealed that Series tractive DC 
motor is very suitable in electric buses transport devices. The starting current was about 460 A, 1s 8 times 
more than the rated current of the motor and 357 A in steady state condition with different loads 
(passengers). In general the Series DC Motor is covering the electric transport requirements and gives the 
best dynamic performance. 
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